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The formation of a two-dimensional electron gas at oxide interfaces as a consequence of polar 
discontinuities has generated an enormous amount of activity due to the variety of interesting 
effects it gives rise to. Here we study under what circumstances similar processes can also take 
place underneath ferroelectric thin films. We use a simple Landau model to demonstrate that in the 
absence of extrinsic screening mechanisms a monodomain phase can be stabilized in ferroelectric 
films by means of an electronic reconstruction. Unlike in the LaAlCU/SrTiOs heterostructure, the 
emergence with thickness of the free charge at the interface is discontinuous. This prediction is 
confirmed by performing first principles simulations of free standing slabs of PbTiCU. The model is 
also used to predict the response of the system to an applied electric field, demonstrating that the 
two-dimensional electron gas can be switched on and off discontinuously and in a non-volatile fashion. 
Furthermore, the reversal of the polarization can be used to switch between a two-dimensional 
electron gas and a two-dimensional hole gas, which should, in principle, have very different transport 
properties. We discuss the possible formation of polarization domains and how such configuration 
competes with the spontaneous accumulation of free charge at the interfaces. 


I. INTRODUCTION 

After the discovery of the formation of a two- 
dimensional electron gas (2DEG) at some oxide 
interfaces^ it was immediately realized that this system 
possessed a number of potential applications. The great 
efforts devoted to the investigation of this unexpected 
phenomenon have indeed yielded a fantastic variety of 
functionalities that can be tailored in these systems such 
as superconductivity^ or enhanced capacitance— Fur¬ 
thermore, the occurrence of a 2DEG in a perovskite sys¬ 
tem opens possibilities for coupling such 2DEG to other 
interesting properties commonly found in different per- 
ovskites, from high-T c superconductivity to multiferroic- 
ity. 

The driving force for the formation of the 2DEG at ox¬ 
ide interfaces, such as the LaAlOa/SrTiOs heterostruc¬ 
ture, is the polar discontinuity at the boundary between 
two materials with different formal polarizations— This 
polarization mismatch has a huge electrostatic cost and 
can favor the formation of free charge that accumulates at 
the interface in order to screen the discontinuity, the pro¬ 
cess sometimes referred to as the “polar catastrophe” — 
One interesting aspect of this phenomenon is that the 
magnitude of the polar discontinuity can be tuned in 
a number of different ways: using different interface 
orientations,- alloying the polar material to effectively 


change its polarization^ or playing with the electrostatic 
boundary conditions of the system .- 9 One possible way 
to manipulate the electrostatic boundary conditions is 
through the coupling with ferroelectricity. This strategy 
has been considered in the pastjiSr— since the sponta¬ 
neous polarization of the ferroelectric material could be 
used to tune the polar mismatch at the interface. In¬ 
deed, the manipulation of the 2DEG in LaAlOa/SrTiOa 
using ferroelectricity has already been achieved experi¬ 
mentally in different ways. In Ref. EH epitaxial strain 
was used to induce a ferroelectric phase transition in 
SrTiOa, whose spontaneous polarization was observed to 
partially screen the polar discontinuity, thus reducing the 
carrier concentration and increasing the critical thickness 
of LaAlO .3 for the formation of the 2DEG. Alternatively, 
in Ref. E2 V. T. Tra and coauthors used a ferroelectric 
over-layer to top-gate the LaAlOa/SrTiOa heterostruc¬ 
ture, being able to induce a metal-insulator transition at 
the interface in a non-volatile way by switching the po¬ 
larization of the ferroelectric. A more radical approach 
is directly to substitute the polar LaA 103 by a ferro¬ 
electric material and use the spontaneous polarization of 
the ferroelectric as the source for the polar discontinu¬ 
ity. This possibility has already been explored from first 
principles.— In Refs. EH and EH it was shown that the 
2DEG could be manipulated with the ferroelectric po¬ 
larization in symmetric KNbOs/ATiOa (A=Sr, Ba, Pb) 
heterostructures. However, the non-stoichiometry of the 
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simulated geometry implied that the interfaces studied in 
those works were metallic by construction. In fact it was 
later demonstrated^ that since centrosymmetric KNI 5 O 3 
is polar with a formal polarization of half a quantum 
of polarization (modulo a quantum of polarization), just 
like LaAlC> 3 , when a [001] interface between this material 
and a non-polar one is grown the ferroelectric polariza¬ 
tion of KNbC >3 tends to compensate the polarity of the 
interface. As a result, the polarization of the KNbC >3 
layer is pinned and its formal value in the ferroelectric 
ground state is approximately zero (up to a quantum of 
polarization) rendering any screening mechanism unnec¬ 
essary. 

Instead, a ferroelectric with a non-polar centrosym¬ 
metric high-temperature phase should be used. In that 
case the polarization in the ferroelectric phase is not 
intrinsically compensated. For a free-standing slab or 
a thin film of such material on top of an insulating, 
non-polar substrate (such as the common SrTiOs), the 
switchable polarization of the ferroelectric could be used 
to manipulate the electrostatic boundary conditions at 
the interface and possibly induce the formation of a 
2DEG. This, of course, would only be possible if such con¬ 
figuration, a monodomain ferroelectric phase screened by 
a 2DEG, is stable, because, unlike LaA 103 , for a ferro¬ 
electric thin film the system has alternative routes avail¬ 
able to minimize or avoid the polarization mismatch. For 
one, in the absence of a screening mechanism other than 
the accumulation of free charge at the interface depo¬ 
larization effects might render the paraelectric configura¬ 
tion as the only stable homogeneous phase of the system. 
But most notably the system can break into polarization 
domains. Strikingly, reports of monodomain phases in 
ferroelectric thin films on insulating substrates are not 
rare in the experimental literatureeven though a 
simple electrostatic analysis reveals that such configura¬ 
tion can only be stable if free charge accumulates at the 
interfaces. In this geometry, interface or surface atomic 
reconstructions, or simple adsorption of ionic species to 
the surface do not provide the necessary screening^ A 
transfer of charge from the surface to the interface or vice 
versa is needed, which, as in the case of LaAlOs/SrTiOs, 
might come from different sources, such as an electronic 
reconstruction or redox processes. In fact, a very re¬ 
cent first-principles study has shown that electronic re¬ 
construction can stabilize a polarization in a BaTiOs thin 
film on top of SrTiOsJ^ In that work, simulations per¬ 
formed for a specific thickness of the ferroelectric film 
showed that, even if the ground state of the system was 
the paraelectric phase, a configuration with a finite po¬ 
larization pointing towards the substrate and a 2DEG at 
the interface was metastable. Neither the polarization 
reversal nor a metal-insulator transition with thickness 
or electric field could be demonstrated, but the work of 
Ref. [lH together with all the previous arguments sug¬ 
gest that the formation and manipulation of a 2DEG at 
ferroelectric interfaces might indeed be possible. 

The complex phenomenology that is expected for these 


systems cannot be explored exclusively within a first 
principles approach. The relative stability of the polar 
configuration with respect to competing phases, the evo¬ 
lution of such competition with the thickness, or the re¬ 
sponse of the system to an external electric field (which is 
the main interest of having a 2DEG in ferroelectric films) 
are issues of prime importance for which a systematic first 
principles analysis is today unfeasible. In this work we 
use a phenomenological model, supported by first prin¬ 
ciples calculations, to confirm that ferroelectricity can 
be used to induce the formation of two-dimensional elec¬ 
tron and holes gases at the interface with non-polar sub¬ 
strates. We discuss the conditions for the stability of such 
configuration, its coupling with external electric fields - 
which gives rise to a discontinuous switching (on and off) 
of the gas-, and its competition with alternative screen¬ 
ing mechanisms such as the formation of polydomain 
phases. The paper is organized as follows: In Sec. [TT] 
we present the model and use it to predict the range of 
stability of a 2 DEG in a prototypical system, in Sec. ITTT1 
we analyze the interaction with an external electric field 
and in Sec. IIVI we discuss the implications of the results, 
in particular how they are affected by the competition 
with polydomain phases and what is the expected phe¬ 
nomenology in the case of the recently proposed hyper- 
ferroelectric materials^ 


II. FORMATION OF A 2DEG AT 
FERROELECTRIC INTERFACES 

A. A simple model 

Phenomenological models have been successfully used 
to rationalize the formation of a 2DEG at polar inter¬ 
faces between paraelectric materials ; 24 : 26 These models 
allow to assess the viability of different processes able 
to screen the polar discontinuity by injecting free charge 
into the interfaces/surfaces. Here we use the same for¬ 
malism to explore the formation of a 2DEG at ferroelec¬ 
tric thin films. We first consider the possible transition 
from paraelectric to ferroelectric with a 2DEG, neglecting 
the competition with the formation of domains (discussed 
in Sec. IIV A[) . 

We start assuming that no extrinsic mechanisms con¬ 
tribute to the screening, thus the electronic reconstruc¬ 
tion (for which electrons from one surface/interface are 
transferred to the opposite one to screen the polarization 
of the film) is the only possible source of free charged 
Throughout this paper, for the sake of conciseness, we 
only speak about 2DEG, but it should be noted that un¬ 
der the assumption of electronic reconstruction the for¬ 
mation of a 2DEG implies the appearance of a corre¬ 
sponding two-dimensional hole gas (2DHG) at the oppo¬ 
site interface or surface. Later on we discuss how this 
model can also be used to describe the basic behavior of 
the system when the free carriers are provided by sur¬ 
face redox processes (like formation of charged defects or 
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FIG. 1. (a) Schematic diagram of the geometry of the sys¬ 
tem. The sign criterion in all the equations throughout the 
paper assume absolute values of both the polarization and the 
free charge density, and their relative orientation is the one 
given in this figure, (b) to (d), schematic band alignment for 
a ferroelectric thin film in various configurations. The corre¬ 
sponding value of the relevant gap, A; the band offset at the 
interface, ij>CB', and the Shottky barrier for electron in the 
presence of an electrode, <j) e , are indicated in each case. 


adsorption of chemical species). 

If the 2DEG forms as a result of electronic reconstruc¬ 
tion, the free energy per unit volume of a ferroelectric 
thin film of thickness d reads 


G=U+^ + ^- + ——— (cr - P) 2 . 


2 gd 2eo 


(1) 


In this expression U is the free energy of the bulk ferro¬ 
electric at zero field, that depends on the polarization P. 
The next two terms on the right-hand side of Eq. [T] ac¬ 
count for the energy cost of promoting electrons from the 
top of the valence band to the bottom of the conduction 
band. The second term of the equation corresponds to 


the cost of the charge transfer across the gap, where a is 
the surface density of free charge, and A is the “relevant 
band gap” of the system in units of a voltage. Neglecting 
surface effects, A is equal to the band gap of the ferro¬ 
electric in the case of a free standing slab [see Fig. QJb)], 
but in general its value depends on the configuration of 
the heterostructure. Fig. |T| illustrates different cases, 
where A is calculated from the band gap of the con¬ 
stituent materials, the band alignment across interfaces 
and the Fermi level of a top electrode (if present). If 
band bending or any other modification of the electronic 
structure occurs at the interfaces, it should also be taken 
into account. Throughout this paper, for the numerical 
estimations we will consider the simplest approximation 
and use the band gap of the bulk ferroelectric for A. The 
third term in Eq. [T] takes into account the energy asso¬ 
ciated to the filling of the bands (the conduction band 
with electrons and the valence band with holes). This 
energy cost is associated with a finite density of states. 
The “reduced density of states’ — used in Eq. [T] is cal¬ 
culated as g = ( fle<7?i)/(<7e + gh) and is expressed in units 
of charge squared per units of area and energy, with g e 
and gh being the density of states (DOS) for electrons 
and holes respectively (we take the densities of states as 
constants, as in a free electron gas in 2D). The last term 
in Eq. |T] corresponds to the electrostatic energy of the 
remnant depolarizing field, and constitutes the driving 
force for the formation of the 2DEG. 

There are several subtleties regarding Eq. [T| that 
should be noted. Firstly, our Landau-type model is re¬ 
stricted to an out-of-plane polarization; it does not in¬ 
clude the possibility of in-plane polarization nor any ex¬ 
plicit strain dependence. Some simple mechanical bound¬ 
ary conditions, such as epitaxial strain, can be implic¬ 
itly accounted for by a renormalization of the coefficients 
in the expansion of [/.— The generalization needed to 
explicitly include these additional degrees of freedom is 
nevertheless trivial: U would depend on the three com¬ 
ponents of the polarization and on the strain, the out- 
of-plane component of P being the only relevant one in 
the depolarization term. Such a model would allow for 
a rotation of the polarization in-plane, which in some 
cases might be a competing mechanism to avoid the elec¬ 
trostatic cost associated to a discontinuity of the out-of- 
plane polarization. Such possibility is however not dis¬ 
cussed here because the in-plane epitaxial compressive 
strain imposed by the substrate hinders the stability of 
an in-plane polarization in the prototypical systems of in¬ 
terest (such as epitaxial PbTiOa or BaTiC >3 on SrTiOa). 
Secondly, in Eq. Q]we assume that the distance between 
electron and hole layers is large enough that exchange 
interactions, as well as excitonic binding, can be disre¬ 
garded. Finally, throughout this paper both a and P are 
the magnitudes of the physical quantities and the signs of 
the different terms in all the equations are valid for the 
geometry and relative orientation of a and P depicted 
schematically in Fig. [Ha)- For an arbitrary sign of the 
free charge with respect to the polarization one should 
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take into account that the term corresponding to the 
band gap energy should read |tr|A. Furthermore, since 
the orientation of the polarization has a strong influence 
on the band alignment of an interface, the relevant gap 
may, in general, have a different value for opposite orien¬ 
tations of the polarization. This difference, together with 
the choice of a suitable thickness, might be exploited to 
switch on and off the 2DEG, as will be discussed below. 

For prototypical ferroelectric materials at zero electric 
field the relevant free energy can be expressed as a Lan¬ 
dau expansion in terms of a single order parameter re¬ 
sponsible for the ferroelectric phase transition 77 , which 
for these materials consists of a soft mode associated to 
a collective shift of the oxygen cage with respect to the 
cations. Alternatively one can use the polarization asso¬ 
ciated with this mode, P v = ^ Z*g 1 where 0 is the unit 
cell volume of the ferroelectric and Z* is the Born ef¬ 
fective charge associated with the mode 77 . The Landau 
expansion of the free energy per unit volume in terms of 
this polarization is expressed as, 

Uo = \a{T - T C )P 2 + \bP* + 0(1*), (2) 

where Tq is the Curie temperature of the material. For 
materials with a second-order phase transition with the 
temperature, the coefficients of P and larger-order terms 
are positive and the energy expansion of Eq. [2] may be 
truncated at the quartic term. This allows most of the 
analysis that follows to be done in terms of analytical ex¬ 
pressions and provides a direct relationship between the 
constants in the Landau expansion and common physical 
properties such as the spontaneous polarization and sus¬ 
ceptibility. A more general discussion should take into 
account higher order terms and, as in the case of im¬ 
proper ferroelectrics, the coupling of polar modes with 
non polar distortions. However here we restrict ourselves 
to materials which can be described by the expression 
in Eq. (2), since this includes some prototypical systems 
such as BaTiC >3 and PbTiC >3 under compressive epitaxial 
strain^! (again, this is the case if either of these materials 
is grown on a SrTiC >3 substrate). In fact, since the phe¬ 
nomenology described in this work is a consequence of 
the “double-well” shape of the free energy as a function 
of the polarization of ferroelectrics, the behavior derived 
from the model should also be qualitatively valid for fer¬ 
roelectrics with a first-order phase transition provided 
that T <C Tc- 

At a given temperature below the phase transition 
(and leaving temperature aside for the time being), we 
rewrite Eq. [2] as 


U 0 


1 

2 £oXt, 


( XP 1 

\4 -P| 



(3) 


where P$ = [(a/b)(Tc — T)] 1 / 2 is the spontaneous po¬ 
larization in the absence of a depolarizing field, and 
£ oXv = [2a(7c ~ T )] _1 is the contribution of 77 to the 
polarizability around P 5 . Xq would correspond to the 


curvature around the minimum of the double well en¬ 
ergy curve Uq(P v ), typically obtained from first princi¬ 
ples performing a series of frozen phonon calculations at 
zero field. The expression in Eq. [3j however, is only 
valid at zero field, since it does not include the extra 
polarization of the electrons and other phonons in arbi¬ 
trary electrostatic boundary conditions. In general, the 
Landau expansion of the energy would be 


U = 




1 ££ 
4 pi 
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2e 0 Xoo e ’ 


(4) 


where the total polarization of the material is 


P = P v + £ 0 X 00 £ = Pri + P e - (5) 

Note that P^ already includes a contribution from the 
electronic cloud contained in Z *, since this is a dynami¬ 
cal charge that takes into account the deformation of the 
electronic charge density with the amplitude of the polar 
distortion at zero field. Accordingly P e and Xoo are the 
extra polarization and susceptibility due to the presence 
of a finite electric held. P e and x '00 account mainly for 
the polarizability of the electronic cloud, thus we will re¬ 
fer to them as electronic polarization and susceptibility 
throughout the paper; however, strictly speaking, these 
two terms also include the contribution of hard modes of 
the lattice . 28 ’ 29 Using the electrostatic boundary condi¬ 
tions of our problem we can express P e in terms of the 
total polarization, P, as 


P e = £ 0Xoo£ = £ 0Xoo- -=Xoo (a-P), (6) 

£0 

which, in turn, can be written as a function of the zero- 
held polarization as 


P = 


Pq T X°°(7 

£(X) 


(7) 


where £qo = Xoo + 1 is the electronic (or background) con¬ 
tribution to the relative permittivity of the ferroelectric. 
Using Eq. [5] through [3 Eq. |T] transforms into 


G = 


1 

2£qX»? 



1 

2£q£oo 


0 




( 8 ) 


Note that the second term in Eq. [ 8 ] looks very similar 
to the last term in Eq. [t] and consequently it could be 
misunderstood as the energy due to the depolarizing held, 
but in fact it contains that contribution as well as the one 
due to the electronic polarization P e . 

Equation [5] can be used to find the equilibrium polar¬ 
ization and surface/interface free charge in a ferroelectric 
thin him. We will neglect for the moment the influence 
of the DOS assuming that g is relatively large. Since 
the DOS term is inversely proportional to g , it decays 
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rapidly for relatively large, but realistic, values of the 
DOS. Therefore, the approximation g —> oo can be made 
without significantly affecting the qualitative behavior of 
the system and simplifying the analysis that follows. Un¬ 
der this approximation, the two equilibrium conditions 
dG/dP v = dG/dcr = 0 yield the system of equations 

5 -^- 7 ^--^)= 0 ( 9 ) 

^5 too 

7 + 7-^-n,) = o. (10) 

a £(Uoo 

under the constraint of a > 0. In the limit of large film 
thickness the solutions of these two equations are cr = P v 
and P v = {—Ps,0,Ps}; and, using Eq. 0 P = P rj . 
This just means that in thick ferroelectric films, the bulk 
tendency dominates, polarizing accordingly, and the free 
charge at the surfaces just follows, precisely canceling the 
depolarizing field (as in shorted boundary conditions). 
As d decreases from oo, P v diminishes, but cr diminishes 
faster, leaving part of the polarization uncompensated 
(the energy cost of the depolarizing field does not com¬ 
pletely overwhelm the energy cost of transferring charge 
from the valence to the conduction band, A). 

In order to study the thickness dependence of the equi¬ 
librium polarization and free charge density, Eq. [9] and 
m can be combined to obtain the condition 


p3 

V 

p3 

US 


Pq_ 

Ps 


= 0 , 


( 11 ) 


where 


1 = 2 Xv~FT- > 


£qA 

Ps 


( 12 ) 


is a characteristic length scale for a given set of parame¬ 
ters. In the plot of P or a versus d , the bulk spontaneous 
polarization Ps defines the scale for P and cr, and the 
length l defines the scale for film thickness, d. 

The evolution of both P and cr is shown in Fig. [2] as 
a function of d. For this plot we used the parameters 
for PbTiC >3 obtained from first principles^, except for 
A, for which we used the experimental band gap, i.e. 
eA = 3.6 eV. The thick film limit (large d) displays what 
was described before, namely, (i) P tends to the bulk 
value Pg, and ( ii ) the free charge density cr tends to 
screen the polarization. The polarization then diminishes 
for thinner films until a critical thickness 


d c — l 


3^ 


(13) 


where a discontinuous jump in all magnitudes occur, and 
below which P = a = 0. The value of the polarization 
at the critical thickness is 

P c = — (14) 

71 V3 ’ [ ’ 

independent of other parameters. The implications of 
Eq. m and El are quite remarkable, in the sense that 



FIG. 2. Total polarization (black solid line), zero-field contri¬ 
bution (dashed line), electronic contribution (solid blue line) 
and surface free charge (solid red line) versus ferroelectric 
film thickness ( d). Ps is the spontaneous polarization of the 
bulk ferroelectric material, d c is the critical thickness for the 
onset of ferroelectric metastability, and dGS is the thickness 
at which the ferroelectric configuration becomes the ground 
state of the system. Inset: Same plot with Xv = 2 and £oo = 1 
to highlight the convergence of cr towards P (in this situation 
P = Pr,). 


these expressions suggest that this is a rather general be¬ 
havior for ferroelectric materials [at least for those which 
respond to the energy expression of Eq. 0, and that the 
main fingerprints of the phenomenon are determined by 
the bulk properties of the material. 

It is interesting to note that the equilibrium screening 
of P by interfacial free charge is quite effective for any 
thickness. The inset in Fig. [2] assumes a value of Xv = 2 
for the plots of a and P versus d, which was chosen for 
illustrative purposes, but represents a very small value for 
any real material. Indeed, a more realistic value (\v = 2 ? 
for PhTiCU) pushes the a curve right onto the P curve, 
as shown in Fig. [2] 

Fig. [3] shows how the solution in Fig. [2] arises. The 
energy functional of Eq. [5] has two possible sets of solu¬ 
tions. For cr = 0 the energy of the system as a function of 
the polarization is a parabola, as in a dielectric material, 
and the equilibrium solution is P v = P = 0. Instead, 
the energy curve for a ^ 0 [solid lines in Fig. 0a)] has 
extrema given by Eq. 1111 The function f(p) = p 3 — p, 
corresponding to the limit of d —> oo of Eq. El and plot¬ 
ted in Fig. 0b), has roots at —1, 0, and 1. As d is re¬ 
duced, the cubic curve shifts upward, then the upper root 
(corresponding to the polarization of the equilibrium fer¬ 
roelectric configuration) diminishes while the middle one 
[corresponding to the energy bump in the G(P ri ) curves of 
Fig. 0a)] becomes positive. The consequence of the this 
is that an energy barrier, which amplitude decays asymp¬ 
totically as d —> oo, separates the paraelectric and ferro- 
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FIG. 3. (a) Energy per unit volume as a function of polariza¬ 
tion for various thicknesses. The numbers next to the curves 
indicate the values of the thickness in each case. Dashed line 
corresponds to a solution with a = 0 while solid lines are the 
curves with a ^ 0. Curves with a ^ 0 possess an equilibrium 
configuration only for d > d c . (b) Graphical solution of the 
equilibrium condition given by Eq. QT] 


electric configurations for every d > d c and the system 
possesses three stable states: zero polarization (insulat¬ 
ing state) and the two opposite polarization orientations 
(2D metallic state). Furthermore, as can be seen in Fig. 
13a), there is a range of thicknesses d c < d < das ~ 1.4d c 
for which the ground state of the system is the paraelec- 
tric configuration, and the ferroelectric one screened by 
the electronic reconstruction is a local energy minimum. 
At the critical thickness d c the minimum of the cubic 
curve touches the axes, meaning that at lower values of d 
the upper root becomes imaginary and ferroelectric con¬ 
figuration ceases to be stable. 

The physical interpretation of the phenomenology de¬ 
scribed above is clear: the appearance of spontaneous 
polarization in the film requires the screening of the depo¬ 


larizing field. This is accomplished by the accumulation 
of free charge that results in the formation of the cor¬ 
responding 2DEG at the interface. Being the screening 
processes a surface effect and the tendency to polariza¬ 
tion a bulk effect, the latter dominates for thick enough 
films, while the former dominates in thin films. Inter¬ 
estingly the transition from paraelectric to ferroelectric, 
or analogously, from an insulating to a conductive inter¬ 
face, is discontinuous. This contrasts with the case of the 
LaA 103 /SrTiC >3 interface for which an effective model 
like this predicts a continuous transition with a gradual 
decrease of the interface charge as the LaAlC >3 thickness 
is reduced. 2 - This of course cannot be observed experi¬ 
mentally since the thickness of LaAlC >3 can only be varied 
in units of the out-of-plane lattice constant, c. The con¬ 
tinuous metal insulator transition in SrTiOs/LaAlOs can 
be shown using an external field to deplete charge from 
the 2 DEGi & 31 i 32 For the ferroelectricity-induced 2DEG 
we will confirm the discontinuous transition under the 
application of an external electric field in Sec. uni 


1. Estimates 


In addition to insights into the character of the tran¬ 
sition, the model, still within the g —> oo approximation, 
allows estimations of the relevant magnitudes. As stated 
above the jump in polarization, Pfj, is Ps/V 3 ~ 0.6 Ps 
and is thus quite universally defined, just dependent on 
the equilibrium polarization of the bulk ferroelectric ma¬ 
terial. 

For the critical film thickness, d c , we can get es¬ 
timates by comparing the results obtained for the 
LaAlOs/SrTiOs interface^! and for the ferroelectric film 
fEos. m\ and fTUll . 


= *^P*oA STO ; = 


3 %/ 3 X ' 

pF 

As 


■e 0 A F , (15) 


where F stands for ferroelectric, and P^ AO refers to half a 
quantum of polarization. Assuming similar values of the 
band gap, and considering now that the critical thickness 
for LaAlOs on SrTiOs is around 4 perovskite layers, and 
that 3V3 ~ 5, 


d F c ~ 20 (x F /x LAO ) layers, (16) 

which can grow quite thick depending on how close the 
temperature is to the bulk ferroelectric T c . 


2. Effect of a finite density of states 

For simplicity, we have assumed so far a large value of 
the reduced DOS, g. This approximation allowed us to 
neglect the energy cost of the filling of the valence and 
conduction bands, simplifying the analysis. The effect of 
a finite DOS is to penalize the accumulation of free charge 
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FIG. 4. Total polarization (black) and free charge (red) ver¬ 
sus thickness calculated using different values for the DOS, g. 
The different approximations for g are g —> oo (solid lines), 
the value corresponding to PbTiOa obtained from its bulk 
band structure^ (dashed), and the two extrema of an esti¬ 
mated range of values for this family of materials (top limit in 
dashed-dotted and bottom limit in dotted lines, respectively). 


and consequently it is expected to shift to larger thick¬ 
nesses the metal-insulator transition. To analyze how 
important the influence of a finite DOS is we plot in Fig. 
[]]the curves of polarization as a function of thickness for 
various values for the DOS, g. For PbTiOa we estimated 
the corresponding 2-dimensional DOS for electrons and 
holes from the bulk band structure obtained from first 
principles. In addition to the usual considerations to cal¬ 
culate a DOS, in these systems one should also take into 
account the fact that the interface lifts some degenera¬ 
cies. We use the results of Ref. HI and [13, which showed 
that after the electronic reconstruction the bottom of the 
conduction band in perovskite titanates has a d xy charac¬ 
ter, to estimate a DOS for electrons of g e /& 2 = 1-2 ■ 10 3 ' 
m _2 J _1 (1.9-10 14 cm _2 eV _1 in more conventional units). 
We performed a similar analysis for the valence band^ to 
obtain the corresponding DOS for holes, which amounts 
to gh/e 2 = 2.5-10 37 m _2 J _1 (4.0-10 14 cm” 2 eV _1 ). Using 
the estimated DOS for PbTiOa we observe that there is 
a significant increase of the transition thickness to about 
~ 1.5 d c , where d c is given by Eq. [13] in the limit of 
g —> oo. The values of d c and dcs predicted by the 
model for PbTiOa, both in the limit of g —> oo and for 
a finite and realistic DOS, can be found in Table HI Ac¬ 
cording to the values in Table [T] the thickness necessary to 
produce a 2DEG in a ferroelectric thin film is large com¬ 
pared with the case of the LaAlOa/SrTiOa interface. For 
ferroelectrics grown under a high epitaxial strain, mech¬ 
anisms for strain relaxation might start playing a role 
in this range of thicknesses, something that would affect 
the parameters in the model. This is not the case for 


TABLE I. Critical thickness for the stability of a 2DEG in¬ 
duced by ferroelectricity, d c , and for the polar phase to be¬ 
come the ground state with respect to a paraelectric configu¬ 
ration, dGS- The thicknesses are calculated both in the limit 
of infinite DOS and using a realistic value for bulk PbTiOar^. 



d c (nm) 

d G s (nm) 

without DOS 

5.5 

7.7 

with DOS 

8.3 

10.5 


PbTiCU on SrTiCU, the system chosen for the numeri¬ 
cal estimations throughout this article, for which epitaxy 
over thicknesses of several hundreds of unit cells can be 
achieved^ Furthermore, the critical thicknesses in Table 
U can be be potentially reduced by an appropriate choice 
of materials, and the model presented here provides a 
simple tool for the screening of optimal systems. 

Taking into account how large the effect of the DOS 
can be, it is worthwhile to explore a range of sensible 
values, since different substrates, terminations or sources 
of free charge (see next Section) could give rise to very 
different values of g. If the source of free charge is, as 
assumed so far, an electronic reconstruction and the ma¬ 
terials involved in the formation of the interface are ox¬ 
ide perovskites, from the estimated values g e and gh for 
PbTiOa we can consider that a given material of this fam¬ 
ily probably displays a DOS in the range between 5 • 10 36 
to 5 • 10 37 m -2 J _1 (8.0 ■ 10 13 to 8.0 • 10 14 cm- 2 eV“ 1 ). 
These two values are used to calculate the two extra 
curves in Fig. [I] Inspecting Fig. 0] one can see that 
in fact materials with a large, but still reasonable, DOS 
might show a dependence of the polarization with respect 
to the thickness very close to the ideal situation. In ma¬ 
terials with a low DOS, on the other hand, the transition 
may take place at thicknesses as large as 2.5 d c — 3 d c . 

The strong influence of a finite DOS on the critical 
thickness for the onset of a 2DEG in ferroelectric films is 
specially noteworthy because such dependence does not 
appear in the case of the LaAlOs/SrTiOs interface. In 
fact, it was demonstrated in Ref. [24] that a model for 
polar interfaces between dielectric materials analogous 
to the one presented here yields a critical thickness for 
the formation of the 2DEG that is independent of the 
DOS (see Eq. 18 in the aforementioned article). 


3. Other sources of free charge 

Although the electronic reconstruction was the first 
mechanism for the formation of the 2DEG in polar inter¬ 
face to be proposed^ alternative processes can also yield 
free carriers to screen the polarization discontinuity.- 24 
Most notably, surface electrochemical processes such as 
adsorption of chemical species or formation of defects in¬ 
duced by the internal electric field in the film are believed 





















to play a major role in the screening at ferroelectric thin 
— and polar interfaces! 26 ’ 42- — To account for this 
kind of mechanism Eq. Q] may be rewritten as 

G = u + 9 j + ^ + k {nQ - F " ? ' (17) 

where C is the formation energy of an isolated redox de¬ 
fect in the absence of an electric field and Q is the charge 
provided by the defect. The defect-defect interaction is 
accounted for in a mean-field approximation by the term 
an 2 /2d. It is easy to show that Eq. [T7] and Eq. |T| are 
in fact equivalent by simply making a = nQ, A = C/Q 
and a = Q 2 / g. Indeed, the analysis presented above for 
electronic reconstruction is parallel to any other planar 
charge screening mechanism associated with an energy 
cost per surface/interface charge. This allows us to treat 
any equivalent screening mechanism with the same equa¬ 
tions, by considering the surface density of free charge, er, 
and taking into account the “effective” gap and DOS rel¬ 
evant for each process. An analysis of the behavior under 
the simultaneous presence of more than one mechanism 
can also be done, as in Ref— for LaAlOa/SrTiOa 


4- Assuming given a or P 


The analysis in the previous subsections assumes equi¬ 
librium, and thus neglects any kinetic effects, which can 
be very important, e.g. in the formation of defects or the 
Zener tunneling of the carriers across the Him. In some 
occasions such kinetic effects may dominate. We can eas¬ 
ily consider the situation in which a certain concentra¬ 
tion of redox defects n has been generated on the sur¬ 
face, e.g. at growth, which are then frozen in. Such sce¬ 
nario corresponds to a ferroelectric layer in open bound¬ 
ary conditions with fixed electric displacement D, where 
D = a = nQ. In this situation a is not a variable but the 
parameter determining the electrostatic boundary con¬ 
ditions of the system. The relevant free energy is now 
simply^ 

G = [/+ J_( ct _P)2. (18) 

2eo 

Using again the transformations for P given by Eq. [6] 
and [7] we get 


G = 



(v-Pn) • (19) 


We can then ask what would be the equilibrium polariza¬ 
tion P for given values of er. Minimizing G with respect 
to P gives the equation 


A 1 f 2 Xy _ -A Pj_ _ 2 Xy a _ q 
-Pg \ £ oo ) Ps £ oo Ps 


( 20 ) 


where the last term is a constant. For any Xy/ £ oo >1/2 
this equation has a single real root, which is positive. 


Considering for simplicity a large value of Xy/ £ cm and 
for values of er not much larger than Pg, the solution can 
be approximated by 


p y 


2 Xy 


1 - 



( 21 ) 


which using Eq. [7] transforms into 


P 


2 w 



( 22 ) 


That is, the polarization responds by compensating 
the effective polarization given by the 2DEG carriers 
(first term) except for a small deviation, which is positive 
(P > cr) for cr < Pg and negative for t 7 > Pg, or, in other 
words, the polarization tends to screen the field generated 
by the fixed surface/interface charge, but with a slight 
tendency towards Pg. Note that under fixed D bound¬ 
ary conditions the free energy of the system scales with 
the volume (it does not have surface terms) and thus the 
behavior obtained is independent of film thickness. Sev¬ 
eral work o 39 ’ 41 have demonstrated that the manipulation 
the surface chemistry can be used to switch the polariza¬ 
tion of a ferroelectric. Furthermore, the first principles 
simulations presented in Ref. |dd] showed how the polar¬ 
ization of a BaTiCU film followed the charge density set 
by charged defects at the surface, as predicted by Eq. [22] 
Similarly one could ask what would be the equilibrium 
concentration of free charge if the polarization P had 
been frozen in by some mechanism. In such (unlikely) 
case, we would minimize Eq.[ 8 ]with respect to a for fixed 
P v and then use Eq. [TJ obtaining 


a 



(23) 


which gives a phenomenology very similar to the cases of 
the LaAlOa/SrTiOa interface, in which the polarization 
is also fixed (in that case to half a quantum) 24 . 


B. First principles simulations 

In order to test the validity of the model, its predic¬ 
tions can be compared with results obtained from first 
principles simulations of ferroelectric thin films. Bear¬ 
ing in mind that first principles simulations are typically 
performed at zero temperature, the results obtained with 
this method should be compared with the low tempera¬ 
ture limit of the model. Nevertheless, as long as the 
temperature is relatively far from the transition one, the 
phenomenology should be qualitatively the same. 

The study presented in Ref. [TU, showing that elec¬ 
tronic reconstruction can stabilize a spontaneous polar¬ 
ization in symmetrical BaTiOs/SrTiCU heterostructures, 
constitutes the first argument supporting the model. In 
addition to this, here we perform additional DFT calcu¬ 
lations on a model system consisting of a slab of PbTiC >3 
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TABLE II. Evolution of the polarization and the energy with 
respect to the non polar phase for PbTiC >3 slabs in vacuum. 
Here we report the values of the polarization of stable struc¬ 
tures after geometry optimizations initialized in a polar con¬ 
figuration, for d < 12 the system spontaneously goes back 
to the paraelectric phase during the relaxation. Energies are 
given per formula unit. 


d (unit cells) 

d (nm) 

P/Ps 

GfE Gpara (m6V) 

10 

4.0 

0 

- 

12 

4.8 

0 

- 

14 

5.6 

0.52 

43 

16 

6.4 

0.62 

32 

18 

7.2 

0.79 

25 


(PbO terminated on both sides) in vacuum, An in-plane 
lattice constant of 3.874 A was chosen to mimic the strain 
of a SrTiC >3 substrate (which is not explicitly included in 
the calculation). Even though such geometry is not rep¬ 
resentative of typical experimental devices, and proper¬ 
ties of the 2DEG such as its confinement or the mobility 
of the charge carriers would be very different from more 
realistic samples like those depicted in Fig. [TJc) and[]Jd), 
we choose here this simple test-case with the sole purpose 
of illustrating some of the basic predictions of the model. 
The calculations were performed within the local density 
approximation, using the siesta code4£ Reciprocal space 
integrations were carried out on a Monkhorst-Pac k 47 i 48 
fc-point grid equivalent to 6 x 6 x 6 in a five-atoms per- 
ovskite unit cell. For real space integrations a uniform 
grid with an equivalent plane-wave cutoff of 600 Ry was 
used. A dipole correction was introduced to avoid spuri¬ 
ous interaction between periodic images of the slab along 
the out of plane direction. Initial coordinates were gen¬ 
erated stacking m unit cells of “bulk-strained” PbTiC >3 
in the ferroelectric phase. Starting from the polar config¬ 
uration we expect that during relaxation the system will 
remain in the metastable configuration predicted by the 
model for d > d c . Then, all the atomic coordinates of 
the slabs were relaxed until the forces were smaller than 
0.04 eV/A. 

In Table El we list the equilibrium polarization and 
energies of the resulting structures as a function of the 
thickness. It was found that for all slabs with d < 12 
unit cells the atoms moved back to the centrosymmet- 
ric positions and the system was insulating. Instead, for 
d > 14 unit cells the slab remained polar and the sur¬ 
faces were metallic. Within the range of thicknesses an¬ 
alyzed here, the energy of the system in the polar phase 
is higher than in the paraelecric one, confirming that the 
ferroelectric/2DEG configuration is still metastable, but 
the energy difference decreases rapidly with increasing 
thickness suggesting that g?gs should be around 10 nm. 

We can now obtain the predicted critical thickness for 
this material according to the model. The four param¬ 


eters, Xv , Ps, A and c, needed to estimate the critical 
thickness were independently computed from DFT cal¬ 
culations on “bulk-strained” PbTiC >3 and found to be 
27, 0.78 C/m 2 , 1.6 V and 4.03 A respectively. The fact 
that the model parameters were obtained from first prin¬ 
ciples calculations means that they also represent the low 
temperature limit and a direct comparison with the first 
principles simulations of the slab can be made. Using 
Eqs. m and □J one gets d c = 2.4 nm ~ 6 unit cells, a 
value significantly smaller than the one estimated from 
the first principles simulations. Nevertheless, if a reason¬ 
able value of the DOS is used (g e /e 2 = 1-2 • 10 37 and 
gh/e 2 = 2.5 • 10 3 ' m _2 J _1 for electrons and holes, as in 
Fig. [1]) the critical thickness increases up to 5.1 nm ~ 13 
unit cells, in excellent agreement with the simulations. 


III. SWITCHING THE 2DEG WITH AN 
EXTERNAL APPLIED FIELD 

The main interest of having a 2DEG in a ferroelectric 
thin Him is that the polar discontinuity at the interface 
might be manipulated in non trivial ways with the appli¬ 
cation of an external electric field. Since the formation of 
a 2DEG relies on the presence of a polarization mismatch 
at the interface, the standard geometry for ferroelectric 
capacitors, with the bottom electrode deposited between 
the substrate and the ferroelectric film cannot be used. 
Instead, in order to preserve the polarity of the interface, 
the electrodes for the manipulation of the 2DEG should 
be placed underneath the dielectric substrate and on top 
of the ferroelectric surface (the latter can be the tip of an 
atomic force microscope), adopting the field-effect tran¬ 
sistor geometry often used to electrically tune the 2DEG 
at the LaAlOASrTiOs interfac e) 31 ! 32 ' 49 (see for instace 
Fig. 1 of Ref. [32|). 

One obvious possible application for this system is the 
non volatile switching of the 2DEG at the interface. If 
the two interfaces or surfaces of the ferroelectric thin film 
are equivalent, switching the direction of the polariza¬ 
tion would simply exchange the 2DEG and 2DHG be¬ 
tween opposite interfaces. If instead, the interfaces are 
dissimilar (if one of them is actually a surface, for in¬ 
stance) switching the polarization might also change the 
effective band gap A, modifying the value of the critical 
thickness d c . If the thickness of the ferroelectric layer 
is close to d c , switching the polarization would then in¬ 
duce a metal-insulator transition at the interface. This 
however requires the application of large electric fields 
to be able to switch the polarization of the ferroelectric 
layer, even larger than for the bulk material since the free 
charge of the 2DEG might respond to the electric field 
and screen it. Nevertheless, as we will see here, switching 
the polarization is not the only way to turn on and off 
the 2DEG at a ferroelectric interface. 

To evaluate the effect of an external electric field in a 
ferroelectric thin film with a 2DEG at one of its inter¬ 
faces we extend the model introduced in previous sec- 
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tions, adding to Eq. [ 8 ] a new term corresponding to the 
interaction of the uncompensated polarization with the 
external field, 

G = u+ ir + ii + k i -’- pf+{ °- p)£ - (24> 

In this expression a positive value of £ represents an elec¬ 
tric field parallel to and with the same sign of the polar¬ 
ization. It should be noted that £ is an external electric 
field and the total field experienced by the ferroelectric 
is <?fe = £ + (a — P)/£ o, with the correct sign criterion. 

The choice of the electric field as the independent vari¬ 
able is in this case natural, since the dependence of the 
polarization on £ only involves the characteristics of the 
ferroelectric layer. In experiments, however, the variable 
that can be directly controlled is typically a gate volt¬ 
age. In such case, the dependence of the polarization 
with the gate voltage requires also knowing details about 
the substrate. Nevertheless, given the specific details of 
a device, the relation between £ and a gate voltage can 
be obtained through 

V = £ (—+d) -^^d, ( 25 ) 

\ £ i / £ o 

where d\ and Ei are the thickness and relative permittiv¬ 
ity of the insulating substrate. With this expression one 
can estimate, for instance, that for a 300 nm thick SrTiCE 
substrate, with a relative permittivity of 300, the maxi¬ 
mum electric field considered in this section (0.4Ps/eo) ws 
requires the application of approximately 30 V between -g 
the top and bottom electrodes. 

... . . . C/5 

The equilibrium polarization of the ferroelectric under 
an applied electric field £ is found after writing Eq. [M] 
in terms of P^, using again Eq. [71 and imposing the 
equilibrium condition dG/dP v = dG/da = 0. To analyze 
the evolution of the polarization and free charge as a 
function of the applied electric field, it is important to 
recall that the energy curves plotted in Fig. [3]are actually 
the result of merging two curves corresponding to two 
different sets of solutions, one for solutions with cr = 0 
and another one for solutions with cr ^ 0. In Fig. [5|a) we 
plot as red solid lines and black dashed lines the curves 
corresponding to cr = 0 and cr ^ 0, respectively, for a 10 
nm thick PbTiCE him and various values of the applied 
electric held. Since Eq. [24] is only valid for cr > 0, the 
sections of the curves corresponding to an electric held 
antiparallel to the screening held due to cr (cr > 0 , £ <0 
or cr < 0 , f > 0 ) are all calculated with cr > 0 and £ < 0 
and reversed with respect to P v = 0 when necessary. 

As discussed before, the energy curve at zero held has 
the appearance of a symmetric triple-well prohle, with 
two metastable ferroelectric states (with 2DEG) and a 
central paraelectric minimum [see curve at the top in 
Fig. [51(a)] . As with the typical ferroelectric double-well 
energy landscape, the effect of the external held is to 
tilt the energy curves, modifying the relative stability of 
the different equilibrium conhgurations and the poten- 


(a) 



I T 1 1 

1 1 ' 1 

D--' 

: 

c 

_ 

A : 

e 

: 

- 

— —‘ ■ _ 

_1_i_1_1_ 

_._1_i_1_ 


_I_i_I_i__i_I_i_I_ 

-0.4 -0.2 0 0.2 0.4 

£ /(P s /e 0 ) 


FIG. 5. (a) Energy curves for a ferroelectric thin film as a 
function of the polarization for different values of the applied 
electric Held (numerical labels next to the curves, in units 
of Ps/eo )■ Red solid (black dashed) sections correspond to 
solutions with cr = 0 (cr ^ 0). Energy curves have been 
shifted vertically for clarity. Light blue dots indicate equi¬ 
librium states as £ is swept. Connected (disconnected) dots 
represent a continuous change (jump) in P v . (b) Hysteresis 
loop for P (black line) and cr (red line) as a function of the 
applied electric field. For the polarization, solid lines repre¬ 
sent the ferroelectric state with the 2DEG while the dashed 
lines correspond to the paraelectric phase. Both (a) and (b) 
are calculated for a PbTiOs thin film with a thickness of 9 nm 
(d > d c ~ 8.3 nm). Points labeled in the top panel correspond 
to those indicated in bottom panel, please find a detailed de¬ 
scription in the text. 
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tial barriers. The energy curves in Fig. Ufa) demon¬ 
strates that the application of an external field can be 
used to switch between the paraelectric (without 2DEG) 
and polarized (with 2DEG) states. Most interestingly, 
the tri-stability of the energy curves in Fig. [5][a) sug¬ 
gests that the metal-insulator transitions should display 
a rather complex hysteresis. In the following lines we 
use the energy curves in Fig. [Sja) to understand the 
shape of the hysteresis loop depicted in Fig. [5{b), ob¬ 
tained using the parameters that correspond to a 9 nm 
thick PbTiOa thin film. As shown in Fig. [5](a), at zero 
electric field both the paraelectric (a = 0, P rj = 0) and 
ferroelectric (cr ~ P v ~ 0.6Ps) configurations are stable, 
with the former being the ground state of the ferroelec¬ 
tric film [point A in panels (a) and (b) of Fig. [5]. As 
we increase the magnitude of the external field, the min¬ 
imum of the cr ^ 0 curve deepens, eventually becoming 
the most stable phase of the system. However there is a 
potential barrier separating the two stable phases, thus 
for small fields (£ < 0.2Ps/e o ), starting in a configura¬ 
tion with a = 0, the system can remain in the paraelec¬ 
tric phase (B). Nevertheless, for a high enough field the 
system eventually overcomes the potential barrier and 
the monodomain configuration as well as the 2DEG are 
switched on (C). If then the electric field is decreased a 
potential barrier prevents the transition back to the para¬ 
electric phase (D). The switching takes place for negative 
fields for which the polar configuration is no longer stable 
(E). For this particular thickness (9 nm), when the en¬ 
ergy curve for P > 0, cr ^ 0 loses its minimum, there still 
exist an energy barrier separating the cr = 0 state from 
the one with P < 0 and cr ^ 0, thus the system remains 
in the insulating phase (E). Therefore, for this thickness, 
the switching of the polarization and surface free charge 
polarity occurs through the non-polar phase. Eventu¬ 
ally, for large enough fields the system switches again to 
a ferroelectric state with a metallic interface (F). 

Assuming the typical geometry depicted in Fig. [He), 
an ideal interface free of defects, and the widely ac¬ 
cepted situation in which free charge at the surface gets 
trapped by defects or adsorbed molecules, the reversal 
of the polarization implies a switching between a 2DEG 
and a 2DHG at the buried interface. This is an inter¬ 
esting result, because a 2DHG has never been observed 
in LaAlOa/SrTiOs. The absence of conductivity at the 
p-type interface is commonly attributed to the fact that 
the polarity of the AlOs/SrO boundary is screened by de¬ 
fects formed during the growth process. For the system 
discussed here, after the deposition the interface is buried 
and protected from further redox reactions. If the ferro¬ 
electric is polarized down as grown, one would expect to 
initially find the 2DEG, which is less susceptible to be 
screened by defects, at the interface, then a switching of 
the polarization might be able to induce the formation 
of the elusive 2DHG. Since electrons and holes can and 
usually do present very different characteristics (such as 
mobility), this possibility constitutes an interesting op¬ 
portunity for the design of new electronic devices based 


(a) 



(b) 




FIG. 6. Hysteresis loop for the polarization (black line) and 
surface free charge (red line) as a function of the applied elec¬ 
tric field for a ferroelectric thin film of two different thick¬ 
nesses: (a) 7.5 nm (d < d c ~ 8.3 nm), and (b) 10 nm (repre¬ 
sentative of the d > d c situation). 


on oxide interfaces, where one could not only play with 
the on and off switching of the 2DEG but also with the 
switching between different gases. 

Another interesting aspect of this system is the fact 
that the shape of the hysteresis loop is strongly depen¬ 
dent on the thickness of the ferroelectric film. Fig. (Ha) 
illustrates the case of a ferroelectric film with a thickness 
below the critical one for the stability of the polar config¬ 
uration at zero field. In this case the switching displays 
two separate loops centered at \£\ > 0. The electric field 
can be used to induce the transition from insulating to 
metallic interface, but the 2DEG would be volatile. This 
hysteresis loop resembles the one corresponding to an 
antiferroelectric, displaying a phase transition to a polar 
state induced by an external electric field and an absence 
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FIG. 7. Phase diagram as a function of thickness d and exter¬ 
nal electric field £ for a PbTi 03 thin film. Circles (crosses) 
indicate an upward (downward) jump in polarization in for¬ 
ward (backward) sweep. Solid (dashed) lines correspond to 
transitions from a ^ 0 (a = 0) to a = 0 (cr ^ 0). Regions 
of paraelectricity and ferroelectricity (with 2DEG or 2DHG) 
coexistence) are shown in dark and light gray respectively. In 
white, different states are accessible depending on the sweep¬ 
ing history. 


of remnant polarization at zero field. Such features max¬ 
imize the electrostatic energy that can be stored in a 
capacitor and suggests that hysteresis loops in ferroelec¬ 
tric thin films might be tailored and optimized for en¬ 
ergy storage applications^ In the opposite limit of large 
thicknesses, shown in Fig. [HDb), once the cr ^ 0 state has 
been reached the switching takes place directly between 
polar states, i.e. between a 2DEG and a 2DHG. 

The whole phase diagram for PbTiC >3 thin films as a 
function of d and £, including the hysteresis regions, is 
shown in Fig. [3 There, the regions of paraelectricity, 
and ferroelectricity (with 2DEG or 2DHG coexistence) 
are shown in dark and light gray, respectively. In white, 
different states are possible depending on the sweeping 
history. For d < d c , at zero field, the ferroelectric state is 
not stable and two separate hysteresis loops are observed 
at finite field, as in Fig. Eta). In Fig. 0 dh marks 
the thickness above which the switching occurs directly 
between the two polar states, like in Fig. Etb). In the 
region between d c and dh all three situations, (cr ^ 0, 
P > 0), (cr 0, P < 0), and (cr = 0), are accessible, 
like in Fig. 0 The phenomenology presented in Fig. 0 
should be quite general, but the shapes of the different 
boundaries between regions of the phase diagram depend 
on the parameters of the material, most notably on the 
DOS, g. 

At this point it is worth noting that the effect of a finite 
DOS in the model cannot be neglected to get the right 
dependence of polarization and free charge with the elec¬ 
tric field. In the limit of g —> oo, starting from the para- 


electric phase and as the electric field increases, the sys¬ 
tem would eventually switch to the ferroelectric/2DEG 
state. However, once in the ferroelectric/2DEG state, if 
the electric field is reversed, since there is no penalty for 
a to grow indefinitely (its dependence with the electric 
field is linear, as in the crA term), the system would never 
switch back to paraelectric or to the opposite polariza¬ 
tion state. In fact, under the g —> oo approximation the 
shape of energy curve of the a ^ 0 configurations does 
not change with the application of an electric field (it 
only shifts up or down), therefore the energy minimum 
corresponding to the ferroelectric state is a metastable 
configuration for any value of the electric field. 

We have assumed here that the free charge can freely 
move from one interface or surface of the ferroelectric to 
the other, implying that no potential barriers are involved 
in these charge transfer processes. This might be reason¬ 
able for very thin ferroelectric films (where tunneling be¬ 
tween the two interfaces is easy) and if the potential bar¬ 
riers for redox reactions at the surface are relatively low. 
Instead, if after the formation of the 2DEG the screening 
charge cannot respond to an external electric field - this 
can be the case of free charge created after an electronic 
reconstruction in a relatively thick film the polarization 
would be pinned by the free charge and the ferroelectric 
would behave as a linear dielectric with a very small sus¬ 
ceptibility. This problem may be overcome by contacting 
the interface with electrodes, puncturing through the fer¬ 
roelectric layer. The metallic contacts could act as the 
source of free charge for the modulation of the 2DEG at 
the interface. 


IV. DISCUSSION 

A. Competition with polarization domains 

In the analysis presented so far we have assumed that 
no polarization domains are formed within the ferroelec¬ 
tric. This competing screening mechanism, however, con¬ 
stitutes the main obstacle for the formation of a 2DEG 
in a ferroelectric thin film. Both the electronic recon¬ 
struction or the redox processes are possible sources of 
screening that can help to stabilize monodomain phases 
in ferroelectric thin films where alternative mechanisms 
(such as metallic electrodes) are not present. Neverthe¬ 
less, the breaking into polarization domains competes 
with the processes discussed here, since in a fully com¬ 
pensated polydomain configuration the net polarization 
charge at surfaces or interfaces is zero, eliminating the 
driving force for an eventual electronic reconstruction or 
surface electrochemical processes. 

Experiments on ferroelectric thin films grown on insu¬ 
lating substrates are abundant in the literature. In the 
following lines we review the most relevant experimen¬ 
tal literature in order to find potential test cases where 
the hypothesis presented here could be investigated. We 
only discuss here those experiments where the ferroelec- 
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trie material is grown directly in contact with the in¬ 
sulator, since the presence of a buffer metallic electrode 
would screen the polar discontinuity at the interface. 


1. PbTi0 3 /SrTi0 3 

The most widely studied ferroelectric heterointerface is 
probably the case of PbTiOa thin films on SrTiCU sub¬ 
strates. In Refs. 0 and 0 , for instance, authors used 
x-ray diffraction methods to determine the polarization 
distribution of PbTiCU films on SrTiC> 3 . In Ref. 0 a 10 
nm film was found to be monodomain as grown with the 
polarization pointing down (i.e. towards the substrate). 
Instead, in Ref. 0 a transition with thickness was ob¬ 
served from polydomain (10 and 20 nm thick films) to 
monodomain with the polarization pointing up (40 nm 
thick films). Refs. 0 and 0 investigate the phase di¬ 
agram of PbTiOs thin films with respect to tempera¬ 
ture and thickness. Ferroelectric films with thicknesses 
ranging from 0.4 (1 unit cell) to 42 nm were grown on 
SrTiCU substrates and the polarization configuration was 
explored using x-ray scattering. At room temperature 
polydomain phases were observed for thicknesses as small 
as 1.2 nm (3 unit cells). For d > 2 nm satellite peaks in 
the x-ray scattering maps disappear but tetragonality is 
consistent with polar PbTiOs, suggesting a transition to 
a monodomain phase. Atomic force microscopy measure¬ 
ments were also consistent with a transition from poly¬ 
domain to monodomain in thicker films.— 

Using coherent Bragg rod analysis D. D. Fong et 
al. were able to obtain a real space mapping of 
the atomic positions (and thus of the polarization) 
in PbTiCU/SrTiCU interfaces.— The authors analyze 
PbTiCU films with thicknesses of 1.6 and 3.6 nm that, 
depending on the cooling process after growth, can be 
stabilized either in a mono (slow cooling down to room 
temperature) or a polydomain phase (fast cooling to 181 
°C). The technique used in this work is sensitive to the lo¬ 
cal electronic density and a priori could be used to image 
the free charge accumulated to screen the polar catastro¬ 
phe at the interfaces. However this would require com¬ 
paring the electronic density of a single sample in the 
monodomain and polydomain configuration (for which 
no electronic transfer is expected). Unfortunately in the 
paper poly and monodomain phases could only be com¬ 
pared in two different samples and any difference in the 
surface termination and interface intermixing obscures 
the possible presence of screening charge. 

Despite the discrepancies about the transition thick¬ 
ness for poly to monodomain phase (that might be as¬ 
cribed to different growth and characterization condi¬ 
tions) the phase diagram that is obtained from this col¬ 
lection of experiments is consistent with the one provided 
in Refs. 0 and0 This phase diagram results from the 
compilation of all the experiments discussed above and 
shows that, at room temperature, for very small thick¬ 
nesses films are paraelectric and as the thickness is in¬ 


creased the ferroelectric films evolve from paraelectric to 
polydomain and then (in a “sluggish transition”, in the 
words of the authors) to a monodomain phase. This sug¬ 
gests that films are polydomain immediately after depo¬ 
sition but became monodomain as they are cooled down 
(unless quenching is used to freeze the poly domain struc¬ 
ture). The driving force for the transition is possibly the 
fact that the orientation of the polarization affects the 
reaction energies of the relevant redox processes at the 
surface, what might make one set of domains more ener¬ 
getically favorable than the other. 


2. BaTi0 3 /SrTi0 3 

In Ref. 0 authors study BaTiCU films grown on 
SrTiCU substrates. SrTiCU exerts a compressive strain 
on BaTi 03 stabilizing a tetragonal phase. In this work, a 
combination of UV Raman spectroscopy and synchrotron 
x-ray scattering is used to test the polarity of films with 
thicknesses of 1.6 to 10 nm. It is found that films of all 
thicknesses are polar at low temperature. At room tem¬ 
perature, films capped with a 10 nm SrTiCU layer were 
polydomain with regular domain sizes. X-ray scattering 
spectra of uncapped films did not show satellite peaks 
but authors were not able to establish whether this was 
due to nonperiodic domains, domain sizes larger than the 
experimental resolution or stabilization of a monodomain 
phase. The model proposed here might offer a possible 
interpretation of this experiment. The uncapped sam¬ 
ples could sustain a monodomain phase screened by the 
accumulation of free charge at the interface and surface. 
In the capped films there would be two possible ways 
for stabilizing a monodomain configuration, (i) Redox 
processes at the SrTiCU surface could provide the nec¬ 
essary free charge to screen the polar discontinuity at 
the ferroclectric/substrate interface, but then the cap¬ 
ping SrTiCU layer should be polarized as well, with the 
corresponding energy cost, (ii) Electronic reconstruction 
within the ferroelectric layer could simultaneously screen 
the polar discontinuity at both interfaces (with the elec¬ 
trode and the capping layer), but the value of A for this 
process is much larger than for the surface electrochem¬ 
ical reactions. In either case, the energy cost of a mon¬ 
odomain polarization in the capped film would be higher 
than for the uncapped one, consistent with the observa¬ 
tion of polydomain phases in the capped samples. 

In Ref. 0 authors report ferroelectricity in 
BaTiCU/SrTiCU/SiCU/Si heterostructures. Even with¬ 
out a top electrode, authors claim that polarization can 
be written in BaTiCU films as thin as 8 nm with the tip of 
an atomic force microscope. For thicknesses of 1.6 nm or 
below, written domains were unstable. They also report 
hysteresis loops in the piezoelectric response without top 
electrode. 

All these experiments, and specially the phase diagram 
for PbTiCU films on SrTiCU substrates provided in Ref. 
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l5ll . consistently support the possibility of stabilizing a 
monodomain phase in ferroelectric thin films grown di¬ 
rectly on insulating substrates. Unfortunately the source 
of screening required to stabilize such configuration is dis¬ 
cussed in very little detail and the possibility of an elec¬ 
tronic reconstruction or alternative process that could 
form a 2DEG or 2DHG at the interface was not contem¬ 
plated in those works. Very recently, however, the for¬ 
mation of a 2DEG at the interface between CaZr 03 and 
SrTiCU has been reported^ Even if CaZrOa is not polar 
in bulk, it is argued in Ref. [52] that the compressive strain 
exerted by the substrate might induce a polarization of 
the CaZr 03 which would be responsible for the forma¬ 
tion of the 2DEG. Although ferroelectricity (switchable 
polarization) could not be demonstrated, the evolution 
of the free charge with thickness is in good agreement 
with the model presented here. Further characterization 
of electromechanical and transport properties in samples 
similar to those discussed in this section, and the com¬ 
parison with the results of the model would be critical to 
assess the hypothesis presented here. 


8.416P 2 

7T 3 £o[1 + (ExEz) 1 / 2 ] ’ 


(27) 


for 180° stripe domains^ The width of domain walls in 
a typical ferroelectric is vanishingly small and remnant 
electric fields in polydomain configurations decay expo¬ 
nentially away from the surface and the domain wall, 
thus, except in the limit of thicknesses of a few unit cells, 
we can approximate |P| ~ Ps throughout the film and 
U = Uq(Ps) as a constant. Using this result and differen¬ 
tiating Eq. [26] with respect to w to find the equilibrium 
domain width for a given thickness of the film, one ob¬ 
tains the well known Kittel law^ 


w 


2 


Ed 

7 


(28) 


Substituting this expression for the equilibrium domain 
width into Eq. [26] leads to the following expression for 
the energy per unit volume 


3. A model for the competition with polarization domains 

The model discussed in Sec. El takes into account the 
balance of the monodomain ferroelectricity with a para- 
electric phase, but cannot say anything about the com¬ 
petition with a polydomain phase. Since experiments 
demonstrate that such a competition exists and the bal¬ 
ance is delicate, here we use a simple model to investigate 
the relative stability of these two phases. 

Since the two screening mechanisms are mutually ex¬ 
clusive we do not need to add new terms to our energy 
expression of Eq. [8] Instead, to analyze the competition 
between the two phases (monodomain and polydomain) 
we need to compare the thickness evolution of the en¬ 
ergy of a thin Him in the two different scenarios. For the 
polydomain phase we assume a 180° domain structure, in 
which straight stripes of the material, all with the same 
width in the direction perpendicular to the domain wall, 
have an out of plane polarization of the same magnitude 
but with alternating orientations. For such idealized ver¬ 
sion of the domain structures typically found in tetrago¬ 
nal ferroelectric thin films the energy per unit of volume 
of the polydomain phase can be expressed as 

G po ly = u + — + Geleo (26) 

where E is the energy per unit of area of a domain wall 
and w is the domain width. The electrostatic energy 
Geiec due to stray fields in the polydomain configuration 
is proportional to the domain width) 53 ' 54 G e i ec = jw/d, 
where the proportionality constant can be calculated to 
be 


Gpoi y (d) = U 0 {Ps) + 2 ■ (29) 

In Eq. [21] the first term is negative and constant, inde¬ 
pendent of the thickness of the film; the second one is 
positive, diverges at d —> 0 (due to the divergence of the 
domain wall density given by the Kittel law) and decays 
with the thickness of the ferroelectric. 

The Kittel law has been shown to be valid down to 
thicknesses of a few nanometers in typical ferroelectric 
thin films^ 6,5 ' This simplified expression for the energy 
is expected to break down below the limit of a few unit 
cells, at which point one should take into account the 
finite width of the domain wall, the stray fields and the 
inhomogeneities of the polarization inside the domains. 
Neglecting these effects, and using again the parameters 
obtained from first principles simulations for PbTi 03 j 2 £ 
we compare the energy of the polydomain configuration 
of Eq. dH (red curve in Fig. [5J with that of a mon¬ 
odomain phase where we allow the possibility of having 
surface charge, as discussed in Sec. Ill A1 given by Eq. [5] 
(black curve in Fig. [5J. As demonstrated before, ferro¬ 
electricity in a monodomain phase becomes stable only 
above a critical thickness, d > d c . In Fig. [8] we are plot¬ 
ting the case where A = C/Q ~ 1 V, corresponding to a 
screening by surface redox processes, that yields a value 
of d c = 1.60 nm. Interestingly, even though for small 
thicknesses the model predicts that the most favorable 
scenario is the formation of polarization domains, it also 
shows that there should be a crossover between the two 
phases, such that above a critical thickness the stabiliza¬ 
tion of a monodomain state by the formation of a 2DEG 
would become energetically favorable over the breaking 
into domains of polarization. For the chosen parameters 
this crossover takes place for a thickness of about 4.6 
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FIG. 8. Energy density of a ferroelectric thin film in mon¬ 
odomain (black curve) and polydomain (red curve) states. 
Curves are calculated for bulk PbTiOa, in the limit of g —► oo, 
using parameters obtained from first principles calculations.— 
In this plot A = 1 V, realistic for surface redox processes. 


TABLE III. Estimated thickness for the mono to polydomain 
crossover, d c _ 0 , in PbTiC >3 thin films. Parameters for the 
model have been obtained from first principles simulations.— 



A (V) 

d c -o (nm) 

Redox processes 

C/Q ~ 1 

4.6 

LDA gap 

1.45 

10 

Exp. gap 

3.6 

64 


nm, well within the range of thicknesses that are typi¬ 
cally grown in ferroelectric thin film experiments. The 
thickness at which the transition from a polydomain to a 
monodomain configuration takes place strongly depends 
on the value of A. This parameter takes different values 
depending on the screening mechanism we are consider¬ 
ing: it is simply the band gap of the ferroelectric if we 
are assuming an electronic reconstruction scenario in a 
free standing slab, and A = C/Q if we consider the pos¬ 
sibility of having electrochemical processes taking place 
at the free surfaces of the material. In the first case it is 
obvious that the band gap underestimation by the tra¬ 
ditional exchange-correlation functionals would greatly 
affect the estimated crossover thickness obtained from 
DFT calculations, as shown in Table uni In the case of a 
ferroelectric thin film on top of a substrate the relevant 
gap A would depend also on the band gap of the sub¬ 
strate material and the band alignment at the interface, 
as indicated in Fig. |T] 

PbTiC >3 was chosen as an example to test the model 
presented here because it is a prototypical ferroelectric 
and has good characteristics to exhibit the formation of 
the 2DEG at reasonable thicknesses, as demonstrated by 
the values in Tables Q] and Mil But in fact one could now 


TABLE IV. Estimated thickness for the mono to polydomain 
crossover in thin films of BiFeOa under a compressive strain 
of 5%. The parameters for the model have been obtained 
from first principles simulations available in the literature (see 
text) except for the relative permittivity, that was calculated 
here = 35, e z = 25). Formation of a 2DEG by means of 
surface redox processes is more favorable than a polydomain 
configurations for all thicknesses. 



A (V) 

dc -o (nm) 

Redox processes 

C/Q ~ 1 

0 

Exp. gap 

3.1 

8.7 


use Eq. [29] and [8] to explore scenarios that could favor 
even further the formation of a 2DEG over the domains 
of polarization. Ideally one would like to find a material 
with large domain wall energy and spontaneous polariza¬ 
tion and small band gap and dielectric constant. A good 
candidate could be the so called “super tetragonal” phase 
of BiFeC >3 (at least from a theoretical point of view, ex¬ 
perimentally the stabilization of these phases requires a 
large in-plane compressive strain and the material usu¬ 
ally forms a mixed phase with other monoclinic phases of 
BiFeOsr^ another disadvantage of this material is that 
experimental samples often present a very large leakage, 
behaving like a semiconductor more than like a true in¬ 
sulator). First principles simulations predict that for in¬ 
plane compressive strains larger than ~ 5% some phases 
of BiFe 03 could display a spontaneous polarizatio n 59 : 60 
of up to 150 C/m 2 and domain wall energies^ of about 
250 mJ/m 2 . A band gap of 3.1 eV has been obtained by 
optical absorption measurements in these highly strained 
phases— All these parameters would yield an estimated 
crossover thickness for the transition from poly to mon¬ 
odomain of 8.7 nm in the electronic reconstruction sce¬ 
nario. Furthermore, if the monodomain polarization is 
stabilized by surface electrochemical reactions this phase 
is more favorable than the polydomain for any thickness, 
as shown in Table IIVI 

This result suggests that even if the breaking up into 
domains is a very effective mechanism for the screening 
of the polar discontinuity at the surface or interface in a 
ferroelectric thin film, the formation of a 2DEG is indeed 
viable and might form for an appropriate combination 
of materials and boundary conditions. The prediction, 
using such a simple model, of a transition from polydo¬ 
main to monodomain as the thickness is increased is in 
good agreement with the experimental observations and 
should constitute a further motivation to explore the sce¬ 
nario proposed here. 

B. Ferroelectric substrate 

If the ferroelectric material is the substrate, and the 
film is a dielectric perovskite the situation is different. 
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FIG. 9. Schematic depiction of the electrostatics involved in 
the case of an insulating thin film on top of a ferroelectric 
substrate, both (a) before and (b) after an electronic recon¬ 
struction. A schematic band diagram of the system is shown 
at the bottom of each panel. Note that for this geometry, 
even before the electronic reconstruction, free charges need 
to be present at the surface to screen the polar discontinuity 
with the vacuum/air region. 


If the substrate is thick and is connected to a metal at 
the back, the polarization discontinuity at the ferroelec¬ 
tric/metal interface would be screened and the ferroelec¬ 
tric would display a strong tendency to develop a finite 
polarization. This, however, gives rise to a polarization 
mismatch at the ferroelectric/dielectric over-layer inter¬ 
face that needs to be compensated to fully screen the 
depolarizing field. In fact, the problem shares many sim¬ 


ilarities with the LaAlOa/SrTiO. 3 . In the latter case, the 
thick substrate imposes a D = 0 electrostatic bound¬ 
ary conditions. Since LaAlOa has a finite polarization 
at zero electric field, the D = 0 condition implies that 
an electric field develops inside the polar material. For 
a thin LaAlOa film, the material polarizes under the ac¬ 
tion of the field, tending to reduce the polar disconti¬ 
nuity. Above a critical thickness though, an electronic 
reconstruction (or a more complex mechanism possibly 
involving redox processes) becomes more energetically 
favorable and a 2DEG forms at the interface. Instead, 
in the case of a non polar dielectric layer on top of a 
thick ferroelectric substrate, the thick ferroelectric im¬ 
poses D = Ps- As in the case of the LaAlOa, in the 
absence of free charge at the interface this condition in¬ 
duces an electric field in the insulating top layer, which 
for small thicknesses will polarize [see Fig. Eta)] un¬ 
til, as the thickness increases, the energy balance favors 
formation of a 2DEG at the interface [see Fig. [Sfb)]. 
There are two important differences with respect to the 
LaAlOa/SrTiOa case though. In the case of the ferroelec¬ 
tric substrate/dielectric film a polar discontinuity exists 
also at the free surface of the dielectric film, which re¬ 
quires the accumulation of some superficial free charge 
(most likely provided by chemical adsorbates, schemat¬ 
ically represented by negative signs in Fig. [9] in accor¬ 
dance with the choice of polarization orientation in the 
ferroelectric) even before the screening at the interface 
sets in. But most importantly, in this system A P can be 
changed (switched) by an electric field and by tempera¬ 
ture. 

Consider for example an SrTiOa film on a BaTiOa sub¬ 
strate. In this case AP is the bulk polarization of the 
substrate, (at low T analogous in magnitude to the polar¬ 
ization mismatch in LaAlOa/SrTiOa), and there would 
be then an instability with a similar critical thickness of 
the film for either electronic reconstruction, the appear¬ 
ance of redox defects, or both. Assuming that beyond 
that critical thickness equilibrium is established in the 
presence of a bulk polarization of BaTiOa parallel to z , a 
2DEG should appear at the interface, of electrons for one 
sign of BaTiOa’s polarization, of holes for the other. If 
switching the substrate ferroelectric, and assuming equi¬ 
libration, the gas of electrons should transform into gas 
of holes and vice-versa. Alternatively one could think of 
a 2DEG (or 2DHG) being switched on and off with T if, 
say, starting from BaTiOa above the ferroelectric critical 
T, and letting it cool down until it polarizes enough to 
give rise to the gas. 

Finally, if one is interested in switching on and off the 
2DEG with a transversal electric field, a polar film can be 
considered, such as LaAlOa, and a ferroelectric substrate 
chosen such that its bulk polarization at the temperature 
of operation is close to the LaAlOa half-quantum. With 
both polarizations aligned, A P ~ 0 and no 2DEG should 
arise. If the ferroelectric is then switched, AP becomes 
approximately a whole quantum, and the 2DEG should 
be strongly populated and stable for quite a thin film, 
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FIG. 10. Total polarization (black solid line), zero field con¬ 
tribution (dashed), electronic contribution (solid blue) and 
surface free charge (solid red) as a function of thickness. d c 
is the one given by Eq. Q3]and in this case it does not mark 
the critical thickness for 2DEG formation. 


possibly allowing for the Zener tunneling to take place 
populating the 2DEG. In this geometry, for either con¬ 
figuration there would still be a polarization discontinuity 
at the free surface, but this would probably be screened 
by redox processes. 


C. Hyperferroelectrics 

Recently a new family of ferroelectric materials has 
been discovered that are predicted to display a finite po¬ 
larization at D = a = 0, hence the reason why they were 
named hyperferroelectrics Their capacity to display a 
finite polarization even in the absence of any source of 
screening might give the impression that the tendency 
of this materials to display 2DEG is even stronger than 
for traditional ferroelectrics. For this reason, in this Sec¬ 
tion we present an analysis for these materials similar to 
the one carried out in Sec. Ill Al Hyperferroelectrics are 
fundamentally equivalent to conventional ferroelectrics in 
the sense that the ferroelectric phase transition is driven 
by an unstable polar mode. Even though their behavior 
in open circuit boundary conditions is radically differ¬ 
ent from that of conventional ferroelectrics, the under¬ 
lying physics is completely analogous. The fundamental 
differences are the low effective charges associated with 
the polar mode and the large polarizability of the elec¬ 
tronic cloud. The consequence of this is that in open 
boundary conditions the electronic polarization screens 
the zero-field contribution and a large fraction of the po¬ 
lar distortion remains stable. As demonstrated in Ref. 
[25l these materials are expected to display a small but 
finite remnant polarization at D = 0. 

Since the basic mechanism for ferroelectricity in hyper¬ 
ferroelectrics is the same as for traditional ferroelectrics, 


the expression of the free energy for these materials is 
again Eq. [5] One can in fact, use that expression to find 
the condition for a ferroelectric to behave as a hyperfer- 
roelectric by calculating the equilibrium polarization for 
o = 0. Using Eq. [TO] with o = 0 and differentiating with 
respect to the polarization, one obtains the equilibrium 
condition 


£lL 1 Eli ( 

PI Ps\ £00 

This equation has solutions 


= 0 


pD —0 
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= 0 ; 


pD —0 
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= ± 



(30) 


(31) 


The non-zero solutions are real only if Soo > 2which 
constitutes the condition for hyperferroelectricity. 

Here we perform an analysis completely analogous to 
the one in Sec. Ill Al to obtain the evolution of polariza¬ 
tion and free charge in a thin film of a hyperferroelectric. 
We use the set of p arameters corresponding to LiBeSb 
reported in Ref. [25| {xy ~ 6 is estimated from the curva¬ 
ture at the minimum of the double well potential). The 
result is plotted in Fig. [TO] Three important differences 
are observed in the curves of Fig. [TO] with respect to the 
corresponding ones for a traditional ferroelectric in Fig. 
[2] First, as expected for these materials, for a thickness 
below the onset for the electronic reconstruction the ma¬ 
terial displays a large polar distortion that is effectively 
screened by the electrons, resulting in a small but finite 
total polarization at D = o = 0. Second, the 2DEG be¬ 
comes stable at a thickness larger than the critical one 
given by Eq. 1131 Finally, neither the total polarization 
nor the free charge display a discontinuous jump at the 
transition. In fact the second and third observations are 
connected. Using Eq. flOl [T3l and 1771 one can demonstrate 
that the value of the free charge at d c is 


Pc 


Ts 

73 



(32) 


This expression yields negative values of 0 for Eoo > 3xy, 
which are not a valid solution of the model. There¬ 
fore for those materials with > 3%^, like LiBeSb, o 
goes to zero (the 2DEG vanishes) in a continuous tran¬ 
sition and for a thickness larger than d c . The range of 
2x v < £00 < %Xy constitutes a third regime in the phase 
diagram with respect to £oo/Xyi where the material is a 
hyperferroelectric but still displays a discontinuous jump 
in the polarization and the free charge. 

As demonstrated in Ref. (2f| the small but finite po¬ 
larization at D = 0 = 0 is a consequence of the strong 
screening provided by the large electronic polarizability, 
but the ferroelectric instability is not necessarily stronger 
than in normal ferroelectrics. We have seen in this section 
that the use of hyerferroelectrics does not favor the for¬ 
mation of a 2DEG as compared with a traditional ferro¬ 
electric, and in fact some of the most remarkable features 
of this system might be lost, such as the discontinuous 
transition of the polarization. 
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V. CONCLUSIONS 


We have used a simple model to demonstrate that 
under the appropriate conditions a monodomain out of 
plane polarization may be stabilized in a ferroelectric thin 
film grown directly on an insulating substrate through 
the formation of a 2DEG at its interface. Although there 
are important analogies with the related polar interfaces 
between dielectric materials, of which LaAlOa/SrTiOa is 
the prototypical example, there are striking differences 
in behavior too. For the 2DEG at ferroelectric interfaces 
the model predicts that a discontinuous transition as a 
function of thickness takes place between the paraelec- 
tric (without 2DEG) and the ferroelectric (with 2DEG) 
phases, with an abrupt jump in both polarization and 
free charge. Also in contrast with the LaAlOa/SrTiOa, 
we have demonstrated that the thickness for this transi¬ 
tion strongly depends on the DOS of the 2DEG. 

One of the key features that was sought in this system 
was the ability to switch on and off the 2DEG as well 
as between a 2DEG and a 2DHG with the application 
of an external electric field. The model shows a com¬ 
plex hysteresis behavior, an effect that poses interesting 


possibilities for energy storage and non-volatile memory 
applications. 

The model has also been used to discuss possible 
strategies to favor this state over the competing para- 
electric and polydomain configurations. We hope that 
the predictions derived from the model, supported by the 
first principles simulations and some of which agree with 
many features from available experiments, will motivate 
the search for 2DEG in these systems. 
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